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HEAT TRANSFER TO BODIES IN A HIGH-SPEED RAREFIED-GAS STREAM 


By Jacrson R. StaLDER, GLEN Goopwin, and Marcus O. Creager! 


SUMMARY 


The equilibrium temperature and heat-transfer coefficients for 
transverse cylinders in a high-speed stream of rarefied gas were 
measured over a range of Knudsen numbers (ratio of molecular- 
mean-free path to cylinder diameter) from 0.025 to 11.8 and 
for Mach numbers from 2.0 to 3.8. The range of free-siream 
Reynolde numbers was from 0.28 to 203. 

The models tested were 0.0010-, 0.0050-, 0.080-, 0.051-, 
0.080-, and 0.12@-inch-diameter cylinders held normal to the 
stream. 

For the case of high-speed flow of a rarefied gas about cylindere, 
the data indicate that fully developed free-molecule flow first 
occurs at Knudsen numbers of approximately 2.0 and the 
temperature-recorery factor depends primarily upon the Knudsen 
number. For Knudsen numbers greater than 0.2, the tempera- 
ture-recovery factor exceeds a value of unity eren though free- 
molecule flow may not be fully developed. 

Over the range of conditions covered in the present tests, the 
Nusselt number ig a function only of the Reynolds number if 
the viscosity and thermal conductivity are based upon stagnation 
temperature and the density based upon free-stream conditions. 
For free-molecule flow, the heat-transfer data are well correlated 
by the theory with the accommodation coefficient equal to 0.90. 


INTRODUCTION 


The subject of heat transfer to or from bodies in a high- 
speed rarefied gas stream has become of importance in 
aeronautics. The recent practical realization of missile flight 
at very high altitudes has stimulated interest in the subject. 
Also, in certain aspects of wind-tunnel testing, notably turbu- 
lence studies in supersonic wind tunnels by means of hot-wire 
anemometers, the effects of slip may not be negligible. (See, 
e. g., reference 1.} 

The various phenomena encountered in high-speed flow of 
gases are often described in terms of the Reynolds number 
and the Mach number of the flow. The former parameter 
may be considered to be a measure of the effect of viscosity 
and the letter a measure of the effect of compressibility on 
the flow field. 

When considering the flow of rarefied gases, however, an 
additional parameter, the Knudsen number, becomes impor- 
tant. The Knudsen number is defined as the ratio of the 
mean-free molecular path to a characteristic body dimension. 
This parameter, which is a measure of the degree of gas 
rarefaction, may also be construed as defining the importance 


1 Supersedes NACA TN 238, “Heat Transfer to Bodies in s High-Speed Rarefied-Gas 
Stream" by Jackson R. Stalder, Glen Goodwin, and Marcus O, Creager, 1951. 


of the effect of molecular motions on the flow field. When 
the Knudsen number is smalt (less than, say, 0.001), the 
effect of the molecular motions on the flow is negligible and, 
in this regime, the gas can be treated as a continuous medium. 
On the other hand, for large values of the Knudsen number 
(of the order of 10), the effect of molecular motions is all- 
important and the phenomena which occur can be completely 
described in terms of the motions of individual molecules— 
in a statistical sense, of course. By virtue of the well- 
established laws of molecular motions, this so-called “free- 
molecule” regime is readily susceptible to analysis and has 
been thoroughly examined (references 2 through 8) with 
respect to both aerodynamic and thermodynamic phenomena. 
Experimental work in this regime with respect. to external 
flows (as opposed to flow in tubes) is, however, quite meager. 
The pioneer work of Epstein in 1924 (reference 3) was 
concerned with the measurement of the drag of spheres; the 
tests were confined to very low-speed motions. Recently, 
the drag and equilibrium temperature of a cylinder in 
high-speed free-molecule flow were reported in reference 8. 

The regime of intermediate Knudsen numbers (0.001 K 
<10) is difficult to deal with analytically and, at this time, 
little success has been attained in attempts to analyze this 
portion of the field of fluid mechanics. Some progress has 
been made, however, by Chapman and Cowling (reference 
9), Schamberg (reference 10), and Grad (reference 11). 

The two general types of approach to the problem used 
by the above authors have been described by Uhlenbeck 
(reference 12). Briefly, the first approach consists of apply- 
ing correction terms to the Navier-Stokes continuum equa- 
tions and boundary conditions, the first approximation 
resulting in a set of higher-order equations—the Burnett 
equations. The other approach is to ¢ttempt to solve the 
Boltzmann equation for the molecular velocity distribution, 


which, in effect, results in obtaining correction terms to the’ 


Maxwell velocity distribution equation. Both of these 
methods are essentially perturbation methods. In the first 
approach solutions are obtained in terms of power series of 
the Knudsen number, while in the second method solutions 
are obtained in terms of the reciprocal of the Knudsen 
number. Both methods will fail for Knudsen numbers of 
the order of unity. f 

As in the case of free-molecule flows, little experimental 
work has been done in the intermediate regime for the case 
of high-speed external flows. However, Kane (reference 
13) has reported results of drag tests on spheres in this 
regime. 
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The general purpose of the research described in this paper 
was to study the heat-transfer processes that otcur between 
a body and a high-speed rarefied gas stream. More par- 
ticularly, it was hoped that the experimental data would 
yield two important results: first, the proper grouping of 
variables to obtain correlation of heat-transfer data in the 
analytically intractable “slip-fiow” regime intermediate 
between the continuum and free-molecule regimes; and 
second, to define approximate values of the parameters 
which delineate the three regimes for the particular con- 
figuration used in these tests. 


Experimenta] data were obtained. on, the heat. trameten 


from circular cylinders in the three regimes of free-molecule, 
slip, and continuum fiows. The data are compared at small 
values of Knudsen number with the data of reference 1 and 
for large values of Knudsen number, with the free-molecule 
analysis of reference 8. 
NOTATION 

Cp specific heat at constant pressure, foot-pounds per 

slug, °F 
d cylinder diameter, feet 
f(a) dimensionless function of s defined by 


f(@)=Z(s'+3)+Z, (2+3) 


g(s) dimensionless function of s defined by g(e)=3 (Z1 +Z) 

g gravitational acceleration, 32.2 feet per second squared 

h heat-transfer coefficient, foot-pounds per square foot, 
second, °F 

I, modified Bessel function of first kind and zero order 

I, modified Bessel function of first kind and first order 


K Knudsen number E ay dimensionless 


k Boltzmann constant, 5.66X10-* foot-pounds per °F 
per molecule F 
l mean-free molecular path, feet 


Af Mach number (=) dimensionless 


Nu Nusselt number ( 42, dimensionless 
p static pressure, pounds per square foot 
Pr Prandtl number (c, £), dimensionless 


O heat transferred from body, foot-pounds per square foot, 
second f 
r temperature-recovery factor, defined by 


T,=T € +> rM 5), dimensionless 
R gas constant, a a per eeu °F absolute 





Re Reynolds number (72 


8 molecular speed ratio (2) an dimensionless 
m 


T  free-stream static temperature, °F absolute 

T. temperature of the cylinder, °F absolute 

T, equilibrium temperature, °F absolute 

U stream mass velocity, feet per second 

vı acoustic specd (yygRT), fect per second 

va most probable molecular speed (V2gRT7), feet per 
second 
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7 mean molecular speed (y SgRI ), feet per second 


Z, dimensionless function of s defined by 


2 
Z:="71 9 (5) exp (-$) 


Z, dimensionless function of s defined by 


OOE 


a ‘accommodation coefficient, dimensionless 

ratio of specific heats, dimensionless 

thermal conductivity, foot-pounds per sccond, square 
foot, °F per foot 

B viscosity, pound-seconds per square foot 


za 


z 


kinematic viscosity (£) square feet per second 


p gas density, slugs per cubie foot 
SUBSCRIPT 
10 refers to the fact that the viscosity and thermal con- 
“ductivity were evaluated at stagnation temperature 
while the density was evaluated at free-stream condi- 
tions 
ANALYSIS 
SOME GENERAL RELATIONS BETWEEN VARIABLES 
Before proceeding with the analysis, it is desirable to set 
down, a few relations which are used throughout. The 
Knudsen number, K=//d, can be related to the Mach and 
Reynolds numbers as follows: ' 


U 
M tbe 5 
hea, Se a a 
x =, 


From reference 14, page 147, the viscosity of a rarefied gas 
composed of hard elastic spheres and having a Maxwellian 
velocity distribution can be expressed in terms of the density, 
mean molecular velocity, and mean-free molecular path, 


p=0.499p01. Thus, very closely, 
M a l2 w l B 
Re Ty K (2) 


The molecular speed ratio s, which is a primary variable in 
the analysis of free-molecule flows, is defined as s= U ftm. 


Since va=V2gRT and v= yygRT, we have a relation — 


between s and M 


E Mo (3) 


HEAT TRANSFER FROM CYLINDERS IN FREE-MOLECULE FLOW 


Heat transfer from a circular cylinder oriented trans- 
versely to a gas stream has been treated, for the case of frec- 
molecule flow, in reference 8. It is useful, however, to 
examine the final equations derived therein with a view to 
putting them in such form as will allow plotting of the 


experimental data on a basis comparable with continuum | 


data. 


HEAT TRANSFER TO BODIES IN A HIGH-SPEED RAREFIED-GAS STREAM 


The equation for heat transfer from a cylinder to a diatomic 
gas stream is shown in reference $ to be (neglecting radiation 
terms) 
2x4? 


DU me 





T gese- Q=0 (4) 


One can define a heat-transfer coefficient per unit area as 


— 2 
hago (5) 





where T, is the equilibrium temperature assumed by the 
cylinder in the absence of heat transfer. The quantity T, 
can be obtained from equation (4) by setting the term Q 
equal-to zero in which case T, equals T, Thus, 


-r fs) 
=P e) È 


Then, substitution of equations (4) and (6) into equation 
(5) yields i 


paPencg(s) (7) 


Oit diffusion 
pumps 


Boundary -layer suction ——-—— 
chamber 






Settling 
chamber 


~~ 
N test chamber 


N poros wall nozzle 


TIT 


If use is made of the relation p/T=geR, equation (7) can be 


rewritten as 





Nua Set repr I@®) 


22°"c, 


Substituting numerical values for the constants, for the — 


case of nitrogen gas, we have, 
l Nu=0.02640RePr ae) de> (9) 


The term g(s)/# is nearly constant for values of ¢ greater than — 


2 and approaches a value of 6yr asymptotically. The 


limiting values are 
[2] ade 
€ -0 & 


Then for values of s>2 we have, very nearly, 


Nu=0.2838aRePr (10) 


For gases, the Prandtl number, Pr, is nearly constant over 
a wide temperature range and, for diatomic gases, has a 









Mechanical vacuum pump 


FictRz 1.—Ames low-density wind tunnel. 
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value of about'0.72 at 0° F. Substitution of this value for 


Pr into equation (10) yields the final equation for heat trans- 
fer from a transverse cylinder in a high-speed (s>2) free- 
molecule flow field 


Nu=0.204al?e (11) 


DESCRIPTION OF EQUIPMENT AND EXPERIMENTAL 
TECHNIQUE 


THE WIND TUNNEL 


The Ames Laboratory low-density wind tunnel, in which 
the present tests were conducted, is an open-jet, nonreturn- 
type tunnel and is described in detail in reference 8. A 
sketch of the major components of the wind tunnel is shown 
in figure 1. 

The nozzles used in the present test, however, differed 
from those described in reference 8 and for that reason will 
be described in detail here. The nozzles used were axially 
symmetric, 1% inches in exit diameter, and had porous walls 
to allow part of the boundary layer to be removed by suction. 
The nozzle contours were determined by calculating the 
shape of the inner core (assumed inviscid) by the method of 
characteristics, as outlined in reference 15, and applying a 
boundary-layer correction to the contour. The boundary- 
layer correction was calculated using the von Kármán 
momentum-integral theory given by Schaaf in reference 16. 
The correction actually applied to the nozzle contours 
corresponded to that calculated for the highest pressure level 
used. 

The nozzles were constructed by stacking 2-inch-diameter, 
thin metal shims to the desired length. The stack of shims 
was held together by four through-bolts and was machined 
to the calculated internal contour. The nozzle walls were 
made porous by discarding every other solid shim and re- 
placing it with small spacers around each of the four through- 
bolts. The shim thickness was varied along the length of 
the nozzle in order to obtain a reasonably uniform suction 
velocity normal to the wall. The nozzle boundary-layer 
removal was accomplished by connecting a chamber around 
the nozzle to one of the main vacuum pumps. 
nozzle is shown diagrammatically in figure 2. 





FIGURE 2.—Cutaway view of porous nozzle. 


The porous | 
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The use of boundary-layer removal resulted in a nozzle 
which produced an essentially constant Mach number over 
the center third of its diameter in contrast with the solid 
nozzles described in‘ reference 8 which produced essentially 
a parabolic Mach. number distribution. Also the use of 
boundary-layer remoral permitted the pressure level of the 
wind tunnel to be changed by a factor of 2%, while the Mach 
number varied only by about +2 percent. <A typical Mach 
number distribution obtained with and without boundary- 
layer suction is shown in figure 3. 


ty © With suction 


o] o Without suction 
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Figure 3.—Varlation of Mach number with distance from center of stream with and without 
boundary-layer suction. Stream static pressure 115 microns of mercury absoluta. 


. The three nozzles used in the tests produced Mach numbers 

of 2.00 4+1% percent, 2.5042 percent, and 3.15+2 percent 
over a range of pressure levels from 80 to 200 microns of 
mercury absolute in the center half-inch of the jet. 


MODELS TESTED 


The models tested were right circular cylinders held, 
normal to the air stream by copper rods fastened to the ends 
of the cylinders. Six models were used in the tests; their 
diameters, test lengths, and types of construction are given 
in the following table: 


Solid, 

Solid, 

Butt-welded thermocouple. 

Solid, 

oe wire wound on coramie 
tube. . 


Fe—Constantan..... 
pian on Pyrex glass- 





Models 1 through 4 had guard heaters wound on the copper 
support rods to allow the temperature of the ends of the 
models to be adjusted in order to prevent end losses. Ther- 
mocouples were soldered to cach end support in order that 
the model end temperature could be measured. Figure 4 
shows the arrangement of the model and its supports in the 
stream. Models 5 and 6 were constructed by wrapping 
0.003-inch-diameter insulated nickel wire on a ceramie tube. 
Voltage leads were soldered to each end of the test length and 
brought out through the ceramic tube to the ends of the 
model. The model was then dipped in a thin solution of 
enamel to obtain a smooth surface. These models were made 
long enough to completely span the nozzle. 
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End 
thermocouples 


Heaters — f 


(electric) 
|| ~~ current 
leads 
“~~ Base 


FIGURE 4.—Dfagrammaile sketch of model and supports. 


Model 3 wes a butt-welded iron-constentan thermocouple 
with the junction located in tbe center of the test length. 

Model 4 was constructed by coating a Pyrex glass rod with 
a very thin layer of platinum. The model was mounted as 
shown in figure 4 by soldering the copper rods directly to the 
platinum coating. 

The temperature of all the models, except number 3, was 
determined from their resistance which was measured by 
comparing the voltage drop across the model to the voltage 
drop across a standard resistor connected in series with the 
model. These voltages were measured with a direct-current 
potentiometer and galvanometer having a least count of 1 
microvolt. The temperature-resistance characteristics of 
each model were determined by calibrating the model in a 
constant temperature bath over the range of temperatures 
encountered in the tests. The temperature of model 3 was 
determined by measuring, with a potentiometer, the thermo- 
couple voltage generated at the junction. 


TEST PROCEDURE 


An impact-pressure tube, a static-pressure tube, and the 
test model were mounted on a common support which could 
he moved through the test area of the jet. The wind tunnel 
was then started and the impact and static pressures were 
measured over the test area of the nozzle. The static- 
pressure tube has been described in reference $ and the im- 
pact-pressure tube was identical with the 0.150-inch-diameter 
impact-pressure tube described in reference 17. Pressure- 
measuring techniques are also described in reference 8. The 
impect-pressure corrections due to viscous effects, given in 


reference 17, were used in order to calculate the stream Mach 


number. 

After the stream conditions had been determined, the test 
model was lowered into the stream and the equilibrium tem- 
peratures of the model test length and ends were measured. 
In the case of models 1 and 2, the guard heaters were used to 
raise the model end temperatures until they were equal to the 
indicated test-length temperature. Heat was not supplied to 
the ends of the other models for the measurement of equi- 
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librium temperature as, for the most part, the equilibrium 


temperatures were, on the average, 10° F below the stagna- 


tion temperature of the stream. 
After the equilibrium temperature of the model was 


measured, the model was heated by passing a direct current 7 


from a battery through the model. For madels 1, 2, and 4, 


the end temperatures were adjusted until there was no 


temperature difference between the ends and the model test 
length. When this point was reached, the voltage across the 


test length of the model and across the standard resistor in_ _ 


series with the model was read on the potentiometer. The 
heat supplied to the model was then the product of the 
voltage drop across the test length of the model and the | 
current through it. 

Model 3, the butt-welded thermocouple, was not used in 
the heat-transfer tests. 


Independent control of the end temperatures was not 


possible with models 5 and 6; however, these models were a 


constructed in such a manner that the resistance to heat. 


flow slong the models was large and the test length was 
confined to the center 3 inch of the models, thus minimizing 
the effect of end losses. 


The radiation heat loss from the models was determined __ 


by measuring the heat loss over a range of model tempera- 
tures, with no flow passing through the tunnel, as the test- 


chamber pressure was decreased toward zero. As it was not — 


possible to obtain a test-chamber pressure below 1 micron of 
mercury absolute, the actual radiation heat loss was obtained 
by extrapolating the measured heat loss to zero pressure. 
Models 1 and 6 were later placed in a chamber in which the 


pressure could be maintained at 0.1 micron of mercury ab- 


solute and the results checked the extrapolated data. 
The resistance of the models at room temperature was 


checked from dey to day and, if a resistance changed enough _ 


so that the temperature of the model was in doubt by more, 


than 2° F, the model was either reealibrated or replaced. 


The tests were conducted at three nominal Mach numbers, _ 
2.0, 2.5, and 3.15, and over a range of pressures from 80 ta ` 


200 microns of mercury absolute. 
The test gas used was commercially pure dry nitrogen. 
A summary of the test data is shown in tabie I. , 
RESULTS AND DISCUSSION 
TEMPERATURE RECOVERY FOR AN INSULATED BODY 


The usual starting point for a set of heat-transfer experi- 


ments at high airspeeds is the determination of the equilib- 


rium temperature of the body; that is, the temperature as- 
sumed by the body in the absence of heat transfer. In the 
case of a flat plate in continuum flow, several analyses are 
available for each of the two cases of laminer and turbulent 
boundary layers. In these analyses, the equilibrium tem- 
perature is expressed in terms of “recovery factor,” r, which 
can be implictly defined as - 


T,=T (1 +71 rit) 


A summary of the available experimental and theoretical 


results is given in reference 18. In brief, the theoretical 
results indicate the recovery factor for laminar boundary 
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« Te obtained from figure 5. . 

è lata taken in University of California tunnel. 
layers to be independent of both Mach and Reynolds num- 
bers and to have a value of Pr’?, All the turbulent-bound- 
ary-layer analyses assumed incompressible flow, hence Mach 
number effects do not appear. The analyses, in general, in- 
dicate that the recovery factor for turbulent boundary layers 
on a flat plate has a value close to Pr'/ with some analyses 
showing a minor dependence of recovery factor on the 
Reynolds number. 

For the case of a cylinder with axis normal to the stream, 
little theory exists and recourse to experiment is required. 
Several experiments have contributed information concern- 
ing the equilibrium temperature of transverse cylinders in 
high-speed air streams, notably those described in references 
1, 19, and 20. There is considerable disagreement among 
the data cited and apparently a definitive experiment with 
a systematic variation of Reynolds and Mach numbers has 
yet to be made. ‘ 

The equilibrium temperature for the case of cylinders in 
free-molecule flow has been determined both theoretically 
and experimentally in reference 8. Good agreement was 
obtained between the experimental results and the somewhat 
surprising theoretical prediction that the cylinder equilib- 
rium temperature would exceed the total or stagnation 
temperature of the gas stream. This phenomenon, of 
course, is in direct contrast to the corresponding phenomenon 
which occurs uņder continuum-flow conditions „where an 
insulated body can have, at most, a temperature equal to 
the stream total temperature and normally does not even 
attain this temperature due to outward heat flow in the 
boundary layer. The anomaly can be explained, however, 
by a@ consideration, in the case of free-molecule flow, of the 
magnitudes of the incident and re-emitted total velocity 
resulting from combination of the stream mass velocity and 
the random thermal velocities. When the total velocity 
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term is squared there results a term, which is twice the 
scalar product of the vector mass velocity and the vector 
thermal velocity. This scalar product term effectively 
increases the apparent internal energy of the gas from the 
continuum value of 3/2 kT per molecule to a value varying 
from 2&T to 5/2 kT depending upon the speed and orientation 
of the body. The apparent internal energy becomes 5/2 kT 
for high speeds and body angles of attack greater than zero. 
This is just equal to the internal plus potential energy per 
molecule of a cube of continuum gas. Therefore, the energy 
incident on the body in free-molecule flow becomes equal to 
that of a continuum for large values of the molecular speed 
ratio for surfaces inclined at angles of attack greater than 
zero. The molecular energy which is re-emitted from the 
surface is'assumed to be equal to the energy of a stream 
issuing effusively from a Maxwellian gas, in equilibrium at 
some yet unspecified temperature, into a perfect vacuum. 
For the case of an insulated body, the temperature of this 
gas is that of the body. The energy of the re-emitted stream 
calculated in this manner is equal tó 2 kT, per molecule. 
The corresponding energy for a continuum gas at the same 
temperature is 5/2 kT, per molecule if both thermal and 
potential energies are considered. For a given re-emitted 
stream temperature, a smaller amount of energy per mole- 
cule is transported from a body for the ease of free-molecule 
flow than is transported in the case of continuum flow. It 
is clear that, if the same total amount of energy, namely 
the incident energy, is to be removed in each case, the 
effusive stream temperature and, hence, the body temperature 


for an insulated body must be higher for free-molecule flow 


than for continuum flow. 4 
The equilibrium temperature for a cylinder in a free- 
molecule flow field has been shown (equation (6)) to be a 


Ratio of measured to theoratical 


Recovery factor, r 
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Frere 3.—Variation of the ratio of measured to theoretical equilibrium tempersture with free-sream Knudsen number for transverse cyHnders. 
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Ficvge 6.—Variation of recovery factor with Knudsen number, 
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function only of the molecular speed ratio, s, or the Mach 
number, Jf. (See equation (3).) 

Because of the different parameters upon which the re- 
covery factor or equilibrium temperature depends in the 
several regimes, it is difficult to plot all the equilibrium- 
temperature or recovery-factor data on a single graph. 
However, for the case of free-molecule flow, the theory pre- 
dicts that the equilibrium-temperature ratio is a function of 
the molecular-speed ratio only. Hence, if we should plot 
the ratio of the measured equilibrium-temperature ratio 
T./T to the theoretical equilibrium-temperature ratio f(s)/g¢(s) 
as a function of Knudsen number, the deviation of the ratio 
from a value of unity should define the lower limit of free- 
molecule flow in terms of the Knudsen number. -These data 
are shown in figure 5 and indicate that the Knudsen number 
must have a value of approximately 2 in order that fully 
developed free-molecule flow exist for a transverse cylinder. 

The recovery-factor data are shown in figure 6, plotted as 
a function of Knudsen number. It can be seen from this 
figure that the data, in the range of Knudsen numbers from 
0.02 to 2.0, are correlated by Knudsen number alone with no 
systematic Mach number effects shown. It is furtber evident 
that, for values of Knudsen number greater than 0.2, the 
recovery factor exceeds a value of unity, presumably due to 
the development of free-molecule-flow effects. Also, it can 
be seen that for values of Knudsen number less than 0.2, 
the recovery factor has a constant value indicating an indc- 
pendence from Reynolds number effects. (See equation (2).) 
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HEAT-TRANSFER-TEST RESULTS 


The presentation of heat-transfer data for these tests is, 


as with the recovery-factor data, beset with difficulty due to. 
lack of a common correlating parameter. The resulta of tho 


heat-transfer tests are shown in one form in figure 7 in which 
Nu is plotted as a function of Hey». The subseript 10 refers 


-to the fact that the viscosity and thermal conductivity wero 


evaluated at tunnel stagnation temperature, while the density 
was evaluated at free-stream conditions. The use of Nuo 
and Rew was suggested by Kovdsznay and Térmarek, who 
found that this particular choice.of conditions for evaluation 
of the gas properties eliminated the need for including ihe 
Mach number as an additional variable. The data of 
Kovasznay and Térmarck (reference 1) are shown as a dotted 
line which was computed from the expression 


Nvy9=0.580 Re? — 0.795 


which is an empirical equation resulting in the best fit of 
their data. All the heat-transfer data obtained during the 
present tests are correlated by the broken line with an 
average deviation of +6 percent. This line was faired 
through the data by the method of least squares and may be 
represented by the relation 


Nuy=0.132 Rey?" (12) 


The solid line represents equation (11) of this paper with the 
accommodation coefficient a having the value of 0.9. This 
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Fioure 7.—Varlation of Nuo with Rew for transverse cylinders in a high-speed stream of rarefled gas 


HEAT TRANSFER TO BODIES IN A HIGH-SPEED RAREFIED-GAS STREAM 


value of æ was chosen to give the best fit of the data obtained 
in the free-molecule range of Knudsen numbers. 

At the other extreme range of Reynolds numbers, how- 
ever, it is seen that the data obtained in the Ames Laboratory 
low-density wind tunnel indicate Nusselt numbers substan- 
tially below those measured by Kovdsznay and Térmarck. 
When these data were obtained, it was suspected that the 
discrepancy was due to blockage effects since the model 
(4-inch diameter) was large compared to the diameter 
(% inch) of that portion of the stream suitable for testing. 
Consequently, some of the tests were repeated using the 
University of California low-density wind tunnel (see refer- 
ence 21 for description) which has a considerably larger 
(approximately 3-inch diameter) testing region. The higher 
Reynolds number data obtained in the University of Cali- 
fornia tunnel again indicated Nusselt numbers below those 
obtained by Kovdsznay and Térmarck although higher than 
those obtained in the Ames tunnel. The fact that the 
present low-density data were obtained in two different wind 
tunnels and (in the case of the Ames test) with three different 
models, while the data of Kovdsznay and Térmarck are 
corroborated by the data of Lowell (reference 20), indicates 
that the difference in Nusselt numbers is real and is not due 
to experimental error.’ 

The discrepancy may be due to one or a combination of 
several effects. First, there is the probability that a differ- 
ence in tunnel turbulence level existed between the tests of 
Kov4sznay and Térmarck and Lowell, which were conducted 


at relatively high pressures, and the present tests which _ 


were conducted at extremely low pressures. It appears 
likely that the turbulence level in a low-density nozzle would 
be low due to the large viscous effects present.’ 

It is also probable that a correction for the effect of vis- 
cosity should be applied to the readings of the static pressure 
tube, as was done in the case of the impact pressure tube 
readings. The value of this correction (which would affect 
the calculated Mach and Reynolds numbers) is unknown 
at present. 

It is apparent that additional tests are required to resolve 
the discrepancy. 

CONCLUSIONS 


The folowing conclusions can be drawn from the data for 
transverse cylinders presented in this paper: 

i. Fully developed free-molecule flow occurs for Knudsen 
numbers of approximately 2 and higher. 

2. For Knudsen numbers greater than 0.2, the tempera- 
ture-recovery factor exceeds unity even though free-molecule 
flow is not fully developed. 

3. For the case of high-speed flow of a rarefied. gas about 
cylinders, the temperature-recovery factor is primarily 
dependent upon Knudsen number. 

4, Over the range of conditions covered in the present 
tests, the Nusselt number is & function only of the Reynolds 


2 It can be observed that the hest-transfer data obtained in the University of California — 


tunnel with the 0.001-inch-diameter model are in good agreement with the data obtained in 
the Ames tunnel; however, it should be pointed out that a shift in the temperature calibration 
of this model resulted in an uncertainty of the order of +10 percent In the data obtained in 
the University of California tunnel. 

3 In passing, it Is noted that McAdams (reference 22) reported large changes in heat-transfer 
rates from cylinders due to vartation of the turbulence level of the alr stream although insuffi- 
cient detalis of the original experimental work were given to enable any conclusions to be 
drawn that would be pertinent to the present discussion. 
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number if the viscosity and thermal conductivity are based 


on stagnation temperature and the density based on free- 
stream conditions. 


5. In free-molecule flow, the heat-transfer data are wel 
correlated by the theory with the accommodation coefñcient 
equal to 0.9. 
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